Two-component mixture of charged particles confined in a channel: melting 
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The melting of a binary system of charged particles confined in a guasi-one-dimensional parabolic 
channel is studied through Monte Carlo simulations. At zero temperature the particles are ordered 
in parallel chains. The melting is anisotropic and different melting temperatures are obtained 
according to the spatial direction, and the different types of particles present in the system. Melting 
is very different for the single-, two- and four-chain configurations. A temperature induced structural 
phase transition is found between two different four chain ordered states which is absent in the mono- 
disperse system. In the mixed regime, where the two types of particles are only slightly different, 
melting is almost isotropic and a thermally induced homogeneous distribution of the distinct types 
of charges is observed. 

PACS numbers; 64.70.dg, 61.46.Hk, 64.70.dj 



I. INTRODUCTION 

Recently, there has been a renewed interest in the 
study of the physical properties of multi-component sys- 
tems [lj43] . The presence of particles with distinct physi- 
cal properties (e.g. size, charge, mass) introduces a com- 
petition between different scales, which is the reason for 
the richer phenomenology in such systems Q. The sim- 
plest multi-component system is a binary mixture of two 
types of particles which, as compared to a mono-disperse 
system [4l4lOj. has been shown to exhibit a very rich 
phase diagram. Such binary mixtures have been studied 
in several different experimental setups such as ion traps 
[ill [l2j , magnetized disks floating at a liquid-air interface 
in an external magnetic field [13,] , dusty plasmas and 
colloidal suspensions [T5| . 

Under specific conditions of density and temperature 
a system of interacting particles will solidify into an or- 
dered arrangement forming a crystal structure. When 
the inter-particle potential is purely repulsive this is of- 
ten referred to as a Wigner crystal, in reference to the 
crystallization of the electron gas predicted by Eugene 
Wigner [l^, and observed experimentally for the first 
time in a 2D system of electrons on the surface of liquid 
helium by Grimes and Adams (T7| . In three-dimensional 
case (3D), a first experimental investi gati on of classical 
plasma crystals can be found in Ref. [l^]- Such an or- 
dered phase appears in several non-electronic systems, 
and has also been experimentally observed, e.g. in ion 
trap s fll| [l^ , colloidal systems [l^ US] : and dusty plas- 
mas |2ll423j [7^ Besides the single-component system, which 
has been widely studied in the last years, crystallization 
may also be observed in multi-component systems of in- 
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teracting particles [H, [3l , and particularly in binary sys- 
tems f?\. Recently, a large number of different equilib- 
rium configurations, which depend on the relative frac- 
tion of the different types of particles, was found in a 
two-dimensional (2D) system of interacting dipoles 6]. 

In the case of a guasi-one-dimensional (QID) system 
of repulsive particles several recent experimental |22h25| , 
analytical and numerical [1, [13, ^B, studies showed inter- 
esting physical properties and promising perspectives for 
practical applications [l^]. Systematic studies of QID 
systems of interac ting particles confined in parabolic @ 
and hardwall [1^ |28| channels revealed many different 
phases of chain-like structures, and even poorly ordered 
configurations. The latter depends on the nature of the 
confinement potential and is favored by hard wall-type 
of confinement. 

An interesting and very important physical phe- 
nomenon is melting. In the 2D case, the transition from 
a solid phase to a complete isotropic liquid phase can 
be characterized by a two-step process with an interme- 
diate hexatic phase [29l - l3]| . However, such a melting 
scenario is not unique and depends on the considered 
system, i.e. the inter-particle interaction, and its dimen- 
sionality. Furthermore, previous papers showed also that 
the melting of 2D clusters is non- universal [H, HI, [11] . 

The specific structural arrangement of the particles in 
confined systems has a very important effect on the na- 
ture and specific process of melting. In the case of 2D 
mono- and bi-disperse confined clusters, structural prop- 
erties of the ordered configuration were found to strongly 
infiuence the melting [sl. [33j. E.g., the melting tempera- 
ture in small clusters of charged particles was found to be 
a nonuniform increasing function of the strength of the 
binding of the cluster [3^]. It was demonstrated in Ref. 
[3^ that the symmetry of the configuration is one of the 
dominant factors that determines the melting tempera- 
ture. In the case of small confined binary clusters, we 
showed recently that thermally induced structural phase 
transitions can occur before melting takes place. In ad- 
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dition, a remarkable temperature induced spatial phase 
separation was found Q . 

In the present paper we use Monte Carlo simulations to 
study the melting of a binary system of classical charged 
particles confined in a gwasi-one-dimensional parabolic 
channel. The structure and normal mode spectrum of 
such a system were recently addressed by us in Ref. lo'] . 
A very rich set of chain-like structures with, in some 
cases, an intrinsic spatial separation of the distinct types 
of particles was predicted. A density dependent band- 
gap in the phonon spectrum was found, and a softening 
of the normal mode frequencies characteristic of a contin- 
uous structural transition from the mixed single-chain to 
the segregated two-chain regime. In the present analysis 
we limit ourselves to an equal density of the two types 
of particles which are assumed to have the same mass 
in order to reduce the number of parameters. Most of 
the results are given for the representative case where 
the charge of one type of particle is twice the charge of 
the other one. Special attention will be given to the par- 
ticular case where the charges are almost equal. In this 
regime at zero temperature, an asymmetrical distribu- 
tion of charges over the chains was found previously [lo| , 
while here we find that temperature leads to a mixing of 
the particles and finally a homogeneous distribution of 
particles. 

The paper is organized as follows. In Sec. II, we de- 
scribe the model and the procedure used to calculate the 
melting temperature. In Sec. Ill the results for the melt- 
ing of different ordered structures corresponding to dif- 
ferent regimes are discussed. Our conclusions are given 
in Sec. IV. 



II. MODEL 

We study a two-dimensional binary cluster consisting 
of an equal number of particles with distinct charges 
Qa and qb, which are allowed to move in the x-y plane. 
The charged particles interact through a repulsive Debye- 
Hiickel (or Yukawa) potential exp{—r/X)/r, are free to 
move in the x direction, and are confined by an external 
one-dimensional parabolic potential which limits the mo- 
tion of the particles in the y direction. A parabolic con- 
finement potential was realized in Ref. [22| for a quasi- 
one-dimensional system of dusty plasma and in Ref. (sBj 
for the case of colloidal particles confined into a droplet 
region. Such a Yukawa inter-particle potential describes, 
e.g. charged colloidal particles in a liquid environment, 
and dusty particles in a plasma. The screening length A 
can be varied by changing the density of the counterions 
(plasma) in the colloidal (dusty plasma) system. In the 
case of a dusty plasma the lateral parabolic confinement 
potential was recently realized by structuring the bot- 
tom electrode into a channel f36j. It was demonstrated 
that the electric field induced potential is parabolic whose 
strength depends on the width of the channel and the 
electric field strength. For the charged colloidal system 



the colloids are e.g. confined between two glass plates 
which can be structured into a channel [2J-|26|. In this 
case the confinement potential can be represented by a 
hard wall potential. It was shown recently [s^l that by 
curving the bottom of the channel it is possible to induce 
a parabolic confinement potential for the colloids where 
the confinement frequency can be varied by changing the 
curvature. 

The potential energy of the system is given by 
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where e is the dielectric constant of the medium the par- 
ticles are moving in, A is the Debye screening length, and 
Tj-^'' = {xi,yi) is the distance of the i^^ particle of type 
A. In order to reveal the important independent param- 
eters of the system, it is convenient to write the energy 
and the distances in units of Eq — (ma;gg^/2e^)^/^ and 
To = (2(7^/mea;g)^/'^, respectively, and define the quan- 
tity a = qb/qa (in the following we express qa and qb in 
units of e, the elementary charge), and the screening pa- 
rameter K = tq/A. In so doing, the expression for the 
potential energy is reduced to 
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and the state of the system is determined now by the 
ratio between the charges a, the number of particles per 
computational unit cell (which is related to the linear 
density), and the dimensionless screening length k. The 
dimensionless linear density n is defined here as the ra- 
tio between the number of particles in the unit cell and 
the length of the unit cell in the unconfined x direction. 
Periodic boundary conditions were considered in the x di- 
rection in order to mimic an infinite system. The largest 
size for the simulation box that we considered along the 
X direction was L = 50 unit cells (N=400 particles), 
and the results were the same as the ones found with 
L 25 (N=200 particles) and L = 12.5 (N=100 parti- 
cles) unit cells, indicating that the results do not depend 
on this length. The temperature is expressed in units of 
To = Eo/kB, where is the Boltzmann constant. From 
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FIG. 1: (Color ) The relative displacements squared (u^) and {Uy) for the distinct types of charges (qa = 1, qt = 2) as a 
function of temperature for different densities (a) n = 0.30, (b) n = 0.50, (c) n = 0.70, and (d) n = 0.74 in the single-chain 
regime. The ground state configuration for the single-chain regime is shown as inset in (b). 



now on, all quantities presented in the manuscript are 
in dimensionless units. The charges, in particular, is in 
units of the elementary charge e. 

The minimum energy configurations and the normal 
modes of this system were systematically analyzed re- 
cently in Ref. [To| . In the present work we employ 
Monte Carlo (MC) simulations to investigate the ther- 
mal dependence of the minimum energy arrangements. 
Starting from a small temperature, typically 5 x lO^'^To, 
the system is equilibrated during typically 10^ MC steps. 
The properties of interest are then measured and aver- 
aged over the next 10^ MC steps. This procedure was 
typically repeated 15-20 times and the quantities of in- 
terest averaged over those 15-20 runs. The temperature 
was increased in steps of the order AT « 10^^ — 10^^, 
and the previous process was repeated. 

It is well-known that in 2D and ID the root mean 
square deviation of the particles from their average lat- 
tice position diverge with the size of the system when the 
temperature is nonzero. This feature was interpreted as 
preventing long-range order in 2D and ID. In the mean- 
time we know that 2D and ID ordering is possible even at 
nonzero temperature. In order to characterize melting, 
we consider therefore a modified Lindemann parameter 
d, m, Lp = {u^)/dl, where (m^) is the difference be- 



tween the mean square displacement of neighboring par- 
ticles from the respective equilibrium positions fg , and dr 
is the mean interparticle distance which is defined below. 
The quantity (u^) is given by the expression 

^ E ]^ E[(^^^ - - (^^'^ - ^^0,)?), (3) 

i=l " j=l 

where () means the average over the MC steps, N is the 
total number of particles in our simulation cell, and Nnt 
is the number of neighbors of the ith particle. Previously, 
it was shown that this modified Lindemann parameter 
characterized the melting transition very well [38]. The 
melting in the x and y directions was characterized by the 
quantities {u'^)/d1 and {Uy)/dy, where the relevant dis- 
tance in each direction, dx and dy , are respectively, given 
by the mean distance between particles in each row, and 
the distance between chains. In the single-chain regime 
we study the behavior of (not Lpy), since there is 
no characteristic distance in the y direction. In the x di- 
rection, the modified Lindemann parameter can 
be defined normally, where dx is the average distance be- 
tween particles along the chain. 

Besides the modified Lindemann parameter described 
above, we also calculated other quantities in order to bet- 
ter characterize the melting phenomenon. An example is 
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the correlation function G,,, given by: 



Gn — 
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dx , dy (defined above) . Tlie 



where rj = x,y and d,j 
function G,, indicates how ordered a configuration is 
Specifically, it compares the distance between particles 
for a given temperature with the distance between par- 
ticles in the ground state configuration (T — 0). In the 
case of a highly ordered configuration, G,, is equal to 
unity. For the completely disordered case the correlation 
function is zero. In the present study, the G^ function 
is calculated for each type of particle separately. This 
means in the single- and two-chain configurations it is 
defined only in the x direction. 



III. RESULTS 
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FIG. 2: (Color ) The slopes jx and 7^ in the low tempera- 
ture interval of the curves (w^) versus T {r) = x,y) obtained 
from simulations in the single-chain regime, and analytically 
as given by Eqs. ([5]) and ((Gjl as a function of density. 



In general, the present binary system can be character- 
ized by an asymmetrical thermal behavior with respect 
to the spatial directions x and j/, and with respect to 
the two types of particles. Anisotropy with respect to 
the spatial direction was also found in the mono-disperse 
system [9j , but we found here that the presence of differ- 
ent types of charges introduces new features. For some 
cases, a very restrictive thermal behavior is observed only 
for one type of particle. 

As a representative example we study the case in which 
the ratio between charges is a = 2, and the screening pa- 
rameter K = 1, which is a typical value for colloids and 
dusty plasmas [s^- The thermal stability of the struc- 
tures will be analyzed for different densities. We will also 
pay special attention to the case when the distinct types 
of charges are almost equal (a ~ 1). Earlier we found 
that in this case, the particles are arranged in a mixed 
configuration where both types of charges are randomly 
distributed over the rows [10[. The zero temperature 
phase diagram of the present system can be found in 
Fig. 3 of Ref. 

We start by presenting the results for the one-chain 
regime, where the mean squared displacements have a 
very particular behavior as a function of temperature. 
Next, the multi-chain regime will be studied, where we 
concentrate on the two- and four-chain regimes. For the 
parameters a = 2, k = 1, and the interval of density 
considered here the system is found only in the mini- 
mum energy structures with one, two, and four chains 
(lOj . Finally, the interesting mixed or disordered case is 
investigated. 

A. The single-chain regime 

The single-chain regime is very peculiar. In princi- 
ple, a Lindemann-like criterion can not be applied here 
since for small densities the difference in the mean square 



displacement behaves linearly with temperature, and no 
clear melting transition is found. For high densities, the 
usual abrupt increase of the mean square displacement is 
found at some critical temperature which is characteris- 
tic for melting, but this is only found along the chain. In 
addition, the dimensionless modified Lindemann param- 
eter presented in Sec. |TT]is not defined in the y direction, 
since there is no characteristic distance in this direction. 
In this case, we still study the thermal behavior of (u^), 
without reference to a melting temperature. 

For a = 2 and n — \ the one-chain configuration [see 
Fig. [TH in Appendix A] is found in the interval of den- 
sity < n < 0.74155. For small densities, n < 0.5, the 
mean square displacements (u^) (and also Lpx) and (w^) 
increase linearly with temperature, typical for harmonic 
oscillations of particles around their respective equilib- 
rium positions [Fig. HJa)]. The slopes of (u^) and {u^) 
are almost the same for the two distinct charges, but 
they are very different along the x and y directions, be- 
ing a factor of almost 20 larger in the former. The reason 
for such an asymmetry can be qualitatively understood 
by the fact that a parabolic confinement in the y direc- 
tion is equivalent to a homogeneous distribution of back- 
ground charge in that direction. This charge distribution 
is constant, independent of temperature, and generates 
a constant pressure in the direction perpendicular to the 
chain of charges. On the other hand, the pressure along 
the chain direction is a consequence of the interaction of 
each charge with all other particles in the chain. Such an 
interaction is strongly affected by the relative positions 
of the particles, i.e. the density, and increases with in- 
creasing density. For small densities the restoring force 
to bring the particles back to their equilibrium position 
is very small which explains the large values obtained for 
{ul). 

For n > 0.5 the behavior of the quantity (Uy) remains 
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FIG. 3: (Color ) The melting temperature of the distinct 
types of charges {qa — 1, qb = 2) and different directions {x, 
y) for the multi-chain regime as a function of density (n). 
The open-square symbols indicate the temperature for which 
a plateau in the (v?) versus T curves appear in the four-chain 
(case 3) regime. 



linear with temperature, but the temperature depen- 
dence of (it^) is qualitatively different. The harmonic os- 
cillation (linear part) observed for small temperatures is 
now followed by a nonlinear regime, which becomes more 
abrupt for larger densities. As illustrated in Figs. [Ijc- 
d), the temperature at which {u'^)/d^ increases abruptly 
becomes smaller with increasing density, and depends on 
the type of particle. For larger densities, fluctuations of 
the particle positions can destabilize the ordered struc- 
ture more easily. This behavior is very different from the 
melting of a classical 2D Wigner solid where the melting 
temperature has the dependence Tc ~ ^/n. 

Another qualitative modification of the (u^) versus T 
curves with increasing density is related with the linear 
regime. Notice that for large densities the slope of the 
linear part of the (m^) curves becomes smaller than the 
one in the y direction. This is a consequence of the fact 
that in the y-direction there is an energy barrier for the 
zig-zag transition , and the barrier height decreases 

with increasing density and consequently the effective po- 
tential minima become shallower. 

An approximately analytic expression for the slope of 
the linear T-dependence of (m^^) can be obtained . Con- 
sider the potential energy given by Eq. ^ , next we limit 
ourselves to the nearest neighbors and then expand the 
Yukawa interaction terms with respect to small oscilla- 
tions. In this way we obtain the deviation of the par- 
ticles, within the harmonic approximation, from which 
we obtain expressions for (u^). Within the harmonic ap- 
proximation we have (u^) — ^ksT from which we found, 
in diniensionless units, the slopes 
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FIG. 4: (Color ) (a) The Lindemann parameter (u^)/d^ (77 — 
x,y) for charges = 1 and gj, = 2), and (b) the correlation 
function Gx as a function of temperature for a system with 
density n — 0.85, a = 2 and k = 1 (two-chain configuration). 
The vertical dotted lines cross the full lines at Gx = 0.5. 



and 



7j, ~ [2 — 2an^(n + k) exp(— k/tt,)] 



(6) 



7a; w [Aan^ K exp{— K / n)] 
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where the density in the single-chain regime is n = 1/a, 
with a the distance between neighbor charges along the 
chain. From Eqs. ([5]) and ([5]) we note that the slope jy 
increases with density. In the x direction, the opposite 
behavior is found for j^- These approximate analytic 
expressions are compared with the results from our sim- 
ulations in Fig. [21 In spite of the rough approximation 
used to obtain Eqs. (O and ([6]), there is good qualitative 
agreement with the numerical results, especially in the 
X direction, which corroborates with the qualitative ex- 
planation commented before. Note that from Eq. ^ we 
find that 7j, diverges for n 0.897 which signals the zig- 
zag transition and compares with the value n = 0.74155 
which we found previously from our MC calculation. The 
difference is due to non-linear effects and contributions 
of particles beyond the nearest neighbor. 



B. The multi-chain regime 

In this subsection the thermal behavior for the two- 
and four-chain structures are presented (Note that a 
three-chain structure is not stable; see Fig. [T^] in Ap- 
pendix A). According to Ref. [l3|, the latter phase can 
be found in two different configurations, defined as case 



1 and case 2. However, we find here that a new four- 
chain (case 3) configuration is possible between these two 
phases. The new phase diagram is shown as Fig. [T2l of 
Appendix A and replaces Fig. 3 of Ref p^ . 

A general feature of the multi-chain regimes is an 
asymmetrical melting in both the direction and the dis- 
tinct types of particles. In the two- and four-chain (case 
2) regimes particles with smaller charge melt first. This 
is not the case for the four-chain (case 1) system. 



1. The two-chain regime 

For n — 0.74155 the system undergoes a continuous 
structural phase transition from the one- to the two-chain 
configuration [l^. For the two-chain configuration, the 
melting is anisotropic with respect to the charges and 
with respect to the spatial directions x and y. The cor- 
responding melting temperatures are summarized in Fig. 
[31 As can be seen in Fig. [3] for both types of parti- 
cles, the melting temperature increases with increasing 
density in the y direction, while the opposite trend is 
found in the unconfined x direction. The explanation for 
such a different dependence is associated with the fact 
that the repulsive Coulomb interaction increases when 
the density increases. Such an increase in the Coulomb 
repulsion leads to a dramatic different behavior in the x 
and y directions. In the latter, and in association with 
the external confinement potential, the amplitude of the 
oscillations are smaller, which explains the behavior of 
the melting temperature in the y direction [Fig. [3]. The 
equilibrium along the x direction is established by the 
balance of the electrostatic force between charges. Ther- 
mal fluctuations of the particle positions destroy such 
balance and, consequently, the ordered arrangement. A 
larger density brings the charges closer to one another, 
increasing the repulsive interaction. In this case, small 
fluctuations of the particle positions (small temperature) 
are enough to destabilize the ordered structure, and to 
melt the system in the unconfined direction. This differ- 
ent behavior along the x and y directions is in line with 
the results found for a single chain as presented in Fig. 

m 

In Fig. [3] we see that for a given density, we find 
in the two-chain regime a smaller melting temperature 
in the y direction. Along the unconfined direction, the 
melting temperature is different for both types of par- 
ticles, being larger for the bigger charge. This result is 
obtained according to the modified Lindemann criterion 
{{u'^)/d^ = 0.1), as iUustrated in Fig. |4Ka) for a system 
with density n — 0.85. This conclusion is corroborated 
by our study through the correlation function Gx defined 
by Eq. (U). As shown in Fig. SKb), similar qualitative 
and quantitative behavior is found for both charges along 
the unconfined direction, where the Gx value for the a- 
particles is slightly smaller which is indicative for a lower 
melting temperature. 

In Fig. [SJ we present the distribution of particles along 
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FIG. 5: (Color ) The distribution of particles along the con- 
fined direction in the two-chain regime (n = 0.85, a = 2 
and K = 1), for different temperatures (a) T = 0.00444, (b) 
T = 0.005, (c) r = 0.006, and (d) T = 0.008. 



the y direction for the same system (density n = 0.85), 
and for different temperatures. Three of those tempera- 
tures, namely T = 0.00444, T = 0.005, and T = 0.006, 
are indicated in Fig. |H[a) by arrows 1, 2, 3, respectively. 
For all temperatures the distribution of particles is asym- 
metric along the y direction, and particles with larger 
charge are more spread out along the confined direction, 
which is due to the larger Coulomb repulsion. For small 
temperature (e.g., T = 0.00444) the distinct types of par- 
ticles are segregated in different chains [Fig. [51(a)], while 
with increasing temperature the two types of particles be- 
come more and more mixed over both chains [Fig. [5jd)]. 
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FIG. 6: (Color) The distances d\/a, d2/a, d^/a, which charac- 
terize the four-chain (case 3) regime, as a function of density 
for a system with k = 1 and a — 2. A sketch of the ground 
state configuration is presented as inset, a is the mean dis- 
tance between particles along each chain. 



The results shown in Fig. [5] indicate that the melting 
in the y direction is not an abrupt process. For T = 0.005 
[arrow 2 in Fig. IH^a)], particles already start to jump be- 
tween chains, which does not occur in the linear regime of 
{ul)/dl [arrow 1 in Fig. IDJa)]. Note that for T = 0.006 
[arrow 3 in Fig. Hl^a)], which is smaller than the melt- 
ing temperature given by the Lindemann-like criterion, 
particles of both chains are already mixed, suggesting 
that the melting process in the y direction is not abrupt, 
and starts with a non-linear increase of the mean square 
displacement. We notice from Fig. IHa) that {Uy^) is 
a rather continuous function while (w^j) exhibits a very 
steep increase at a certain temperature. 

Along the chain direction, the temperature for which a 
deviation from the linear decay of Gx occurs [Fig. IHJb)], 
coincides with the one along the x direction for which the 
Lindemann parameter deviates from a linear temperature 
behavior [Fig. SJa)]. Note that in this case, the devia- 
tion of the Lindemann parameter from the linear regime 
starts at the same temperature (T = 0.007) for both 
types of charges [Fig. HJa)] as observed in Fig. |3Jb) for 
Gx . In the unconfined direction, the correlation function 
together with (u^J suggests a continuous disordering of 
the system with increasing temperature. 



2. The four- chain (case 1 and case 3) regimes 

A further increase in the density brings the system 
to the four-chain regime, which is found in three differ- 
ent configurations. In all of them (case 1, case 2, case 
3), the different types of particles are symmetrically dis- 
tributed with respect to the a:-axis, but segregated in 
distinct chains. In this section we discuss the thermal 
behavior of case 1 and case 3. 
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FIG. 7: (Color) (a) The Lindemann parameter (u^)/d^, and 
(b) the correlation function Gr, [r] = x, y) for charges qa — 1 
and Qf, = 2 as a function of temperature for a system with 
density n = 1.60, a = 2 and k = 1. Particle trajectories 
over 100 MC steps for temperatures T = 0.001188 (arrow 1) 
and T = 0.001875 (arrow 2) are shown as insets in (a) and 
(b), respectively. Black (red) circles represent particles with 
charge qa ~ 1 (qt ~ 2). (c) The distribution of particles in 
the y direction, (d) The pair-correlation function along the 
X direction, a is the mean distance between particles along 
each chain in the Eround state confiEuration. 
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In the four-chain (case 1) regime, which is found in 
the density interval 0.9 < n < 1.16, chains with equal 
charges are displaced with respect to each other by a/4 
(a is the distance between particles in each chain) along 
the chain direction (a;-axis), while neighbor rows with 
distinct types of particles are displaced by a distance 
a/2 along the a;-axis. The distance between the inter- 
nal chains (consisting of particles with the same charge) 
is larger than the distance between the internal chains 
and the external ones [l3|- This is interesting because 
the interaction between chains with distinct charges is 
intuitively expected to be larger than the interaction be- 
tween chains with the same lower charge. In the case 1 
regime the melting temperature in both spatial directions 
and for the two different particles is almost the same as 
shown in Fig. [3] 

For 1.16 < n < 1.68 the system suffers a continuous 
transition to the four-chain (case 3) regime which was 
not predicted in Ref. [l3|. Details of the internal con- 
figuration of particles in such a regime is presented as 
inset in Fig. |51 The main difference with respect to the 
case 1 regime is the density dependence of the distances 
(di, ^2,(^3) which characterizes the ground state configu- 
ration. In the four-chain (case 1) regime di — d2 — a/4, 
and ds = a/2 (a is the distance between particles along 
each chain). As will be presented in the next section, 
the four-chain (case 2) configuration is characterized by 
di = d2 = ds = a/2. 

The behavior of the melting temperature in the four- 
chain (case 3) regime is rather distinct when compared to 
the one in case 1 . There is a sudden increase of the melt- 
ing temperature of the external chains in the y direction 
for n > 1.4. On the other hand, the melting temper- 
ature for the internal chains decreases with increasing 
density. The same decrease in the melting temperature 
is observed for the external chains in the x direction. The 
reason for the modified density dependence of the melting 
temperature is a thermally structural phase transition in 
which the system changes from the four-chain (case 3) 
to the four-chain (case 2) configuration. For example, in 
Fig. El^a) the Lindemann parameter (it^^)/(i^^ (77 = x,y) 
for charges qa and qb is presented as function of temper- 
ature for a system with n — 1.60. As can be seen, there 
is a small plateau in the mean square displacement for 
both charges, but only in the unconfined direction. As 
observed previously in the case of circularly confined fi- 
nite size clusters [3J|, the plateau is a typical signature 
of a thermally induced structural phase transition, and 
we find here that the plateau is, in fact, associated with 
such a transition. In the insets in Figs. El^a) and[71[b), 
we present typical trajectories of the particles obtained 
over 100 MC simulation steps for temperatures below 
and above the plateau [indicated by the arrows 1 and 
2 in Fig. Wi^)]- Around the structural phase transition 
the system is in the harmonic regime, since the mean 
square displacement increases linearly with temperature. 
For temperatures before the plateau, the particles are ar- 
ranged in the four-chain (case 3) regime, while above the 



plateau the trajectories of the particles indicate that the 
system is in the four-chain (case 2) regime [lOj], whose fea- 
tures are described in the next section. Such a thermally 
induced structural phase transition is observed only for 
densities n > 1.50 in the four-chain (case 3) regime. The 
temperature for which the plateaus appears as a function 
of density is indicated in Fig. [3] by the open square sym- 
bols. Above those temperatures the system is found in 
the case 2 regime, as indicated in Fig. |3l 

From the analysis of the Lindemann parameter, we find 
that the thermal behavior depends on the spatial direc- 
tion, which is qualitatively similar for both types of par- 
ticles. However, if the correlation function (77 = x, y) 
defined in Eq. is considered [Fig. EJb)], we still ob- 
serve a similar thermal dependence in the x direction for 
the charges qa and qb- But this is different in the y di- 
rection. The thermal dependence of Gy for charge qa 
(internal chains) suggests a continuous loss of ordering 
with increasing temperature, and differs from the one of 
particles with charge qb (external chains). The latter thus 
indicates an ordered structure due to the balance between 
the Coulomb repulsion from the charges in the internal 
chains and the external pressure from the parabolic con- 
finement potential. 

For high temperatures, e.g. T = 0.0277, the inter- 
nal chain-like structure of the system is disappeared, 
while the external chains still remain well defined [Fig 
[Tljc)]. However, for both internal and external chains 
a liquid-like behavior is found along the chain direction. 
The mean distance between particles with smaller charge 
{qa) in the x direction becomes half the distance be- 
tween particles with larger charge (qb). This is also illus- 
trated by Fig. E^d), where the pair correlation function 
gix) = 'E;=ij{S[x - {x^ - Xj)]) is shown. 

3. The four- chain (case 2) regime 

For T = there is a continuous structural phase tran- 
sition, at n « 1.68, from case 3 to case 2 (see Fig. [S]). 
In the four-chain (case 2) regime, chains with the same 
charge are displaced by a/ 2 with respect to each other 
along the x direction. The density dependence of the dis- 
tance between chains has an opposite behavior to the one 
of case 1, i.e. the distance between chains with distinct 
charges is larger than the distance between the internal 
rows [lol |. The different internal structure of the four- 
chain (case 2) regime leads to a different temperature 
dependence for the mean square displacements. As can 
be observed in Fig. [31 the melting temperature for both 
types of particles, and in the two directions has a differ- 
ent behavior as a function of density in the two phases 
of the four-chain regime. 

The particular arrangement of the charges [see inset in 
Fig. [6] results in a larger distance between the internal 
(qa) and the external (qb) chains than that found in case 
1. As a consequence, the external confinement potential 
acts more strongly on the particles placed in the outer 
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FIG. 8: (Color) (a) The Lindemann parameter {u^)/d^, and 
(b) the correlation function Grj {r] = x,y) for charges qa and 
qb as a function of temperature for a system with density 
n = 2.40, a = 2 and k = 1. (c) The distribution of particles 
in the y direction, (d) The pair-correlation function along the 
X direction, a is the mean distance between particles along 
each chain in the ground state configuration. The peaks of 
the pair correlation function for the b-particles is fitted by an 
algebraic decay represented by the dotted blue curve. The 
legend in (b) is the same of (a), and the legend in (d) is the 
same of (c). 



rows, resulting in a larger melting temperature in the 
y direction. Fig. |3] indicates a clear anisotropy in the 
melting temperature of the external chains with respect 
to the X and y directions. 

Now, considering the inner chains we observe that par- 
ticles in different rows are already close such that fluctua- 
tions of the particle positions are sufficient to destroy the 
ordering in both directions. This fact occurs for a tem- 
perature smaller than those needed to melt the external 
chains. We can conclude that according to the modified 
Lindemann criterion considered in the present paper, the 
melting in the four-chain (case 2) regime is anisotropic 
with respect to both spatial directions and the distinct 
types of charges. 

We resort again to the correlation function G-q {rj = 
X, y) and other quantities in order to better understand 
the temperature dependence. As an example, in Figs. 
|H{a) and IHJb), a comparison between the temperature 
dependence of the Lindemann parameter and the corre- 
lation function G,, for density n — 2 AO is made. In both 
the X and y directions the Lindemann parameter has the 
same qualitative behavior as a function of temperature 
for both types of particles, i.e. a linear behavior fol- 
lowed by an abrupt non-linear increase of {u'^)/d^ [Fig. 
IUJa)]. On the other hand, G,, presents a very different 
temperature dependence in the x and y direction, and 
for both types of charges. As in the case 1 regime, Gy 
shows a continuous loss of ordering as a function of T 
for the internal chains {qa), which is not observed for 
the external ones. In the x direction, the behavior of 
the internal and external chains is similar only for small 
temperatures (T < 0.01). For T > 0.01, the behavior 
of Gx for the internal chains exhibits an abrupt decay 
in the corresponding temperature interval associated to 
the melting according to the Lindemann parameter cri- 
terion. For particles in the external chains, we find a Gy 
T-dependence that suggests a slowly disordering with in- 
creasing temperature, which is not seen from the Linde- 
mann parameter [Fig. [SJ^a)]. For example, at the tem- 
perature T — 0.037 the Lindemann parameter indicates 
that the system is completely melted, but Gy indicates 
that only the internal chains are completely disordered 
(i.e., Gx^Gy^ 0). 

In Figs. [SJc) and Efd), we show the distribution of 
particles along the y direction and the pair correlation 
function along the unconfined direction for T — 0.037, 
respectively. As can be observed, the initial two-chain 
structure of the internal chains is completely destroyed 
in the y direction, while it is relatively preserved in the 
external chains. Along the x direction. Fig. [Sfd) indi- 
cates that the internal chains are completely disordered, 
but the pair-correlation function for particles in the ex- 
ternal chains presents an algebraic decay, which is indica- 
tive of long-range order. Notice that the Gx function for 
the b particles still presents a finite value for T = 0.037 
in agreement with the g{x) result. For very high tem- 
peratures, and as found in the case 3 regime, the mean 
distance between the particles {x direction) with smaller 
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FIG. 9: (Color) The Lindemann parameter (u^) /d^ for a sys- 
tem in the mixed regime {qa = 1 and 55 = 1.02) as a func- 
tion of temperature for different densities (a) n = 2.20, (b) 
n = 2.50, and (c) n = 2.80. 



charge becomes half the one for larger charge (first peak 
in the pair correlation function), which indicates that a 
disordered one chain is found in the internal part of the 
system. The results presented here indicate a complete 
anisotropic (direction and type of charge) melting for the 
binary system. 



C. Mixed regime 

In the almost mono-disperse case (i.e. a k, 1), the 
minimum energy arrangement is the three-chain config- 
uration with the charges randomly distributed over the 
chains. Neighbor chains are displaced by a distance a/2 
with respect to each other along the unconfined direction 
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FIG. 10: (Color) The distribution of particles {qa ~ 1 and 
qt = 1.02) along the confined y direction for a system in 
the mixed regime for temperatures (a) T = 0.005, and (b) 
T = 0.01083. 



[inset of Fig. [TT]]. For k = 1, the interval of density in 
which the mixed regime can be observed depends on the 
ratio between the size of the charges a UM- Here, we 
focus on the representative case a — 1.02, i.e. particles 
have a difference in their charges of 2% where, for T = 0, 
the disordered state was found in the density interval 
2.14<n<3.23. 

The temperature dependence of the Lindemann pa- 
rameter is presented in Fig. [9] for systems with a = 1.02, 
K = 1, and densities n — 2.20, 2.50,2.80. In spite of the 
small density interval in which the mixed regime is found 
[lo| , a different qualitative behavior is observed with in- 
creasing density. In general, the melting temperature 
increases with increasing density. In addition, the melt- 
ing temperature in the y direction becomes distinct for 
both types of particles. In the x direction, the melting 
temperature depends on the density, but it is the same 
for the two set of particles. 

The random arrange in the mixed regime generates an 
unbalanced distribution of the different charges over the 
chains, as can be seen in Fig. [TUl a). where the distribu- 
tion of particles along the confined direction is shown for 
the system with density n = 2.20 {a = 1.02, k = 1) at 
T — 0.005. An interesting feature is observed when the 
temperature is increased. In this case, a more homoge- 
neous distribution of both types of particles is induced, 
where particles with larger charge are located predomi- 
nantly at the edge chains [Fig. [TUf b)]. A similar behavior 
is found for all densities in the mixed regime. It is inter- 
esting to comment on the more homogeneous distribution 
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FIG. 11: (Color) The melting temperature of the distinct 
types of particles {qa = 1 and qi, — 1.02) as a function of 
density for a system in the mixed regime. An example of the 
ground state configuration is shown as inset 



of particles which can also be noticed from the correla- 
tion function. In particular, we present as inset in Fig. 
inia) the correlation function Gx as a function of temper- 
ature for the system with n ~ 2.20 {a — 1.02, k — 1). 
For T > 0.005 the correlation function Gx is larger for 
particles with larger charge {qb — 1-02), which coincides 
with the rapid increase of {Uy) at which particles start 
to jump between chains leading to a more homogeneous 
distribution of the different particles over the chains. At 
this point the system melts. The density dependence 
of the different melting temperatures of the system with 
a = 1.02 and k = 1 in the mixed regime is summarized 
in Fig. HH 

IV. CONCLUSION 

We studied melting and thermal induced transitions 
of a two-dimensional binary system consisting of an 
equal number of particles with different charges interact- 
ing through a screened Coulomb potential, and confined 
along one direction by a parabolic potential. The differ- 
ent phases (see Fig. [T^ that nucleate as a function of 
density jTo| exhibit a very distinct and nontrivial tem- 
perature dependence. 

The melting temperature was determined from a mod- 
ified Lindemann criterion Q. Other quantities as the 
pair-correlation function, the distribution of particles, 
and the translational correlation function were in addi- 
tion studied to corroborate the obtained values for the 
melting temperature and to enhance our understanding 
of the melting process of the different ordered structures. 

Wc found that melting is anisotropic with respect to 
the spatial x and y directions, and with respect to the 
distinct type of particles. In the multi-chain regime, de- 
pending on the internal structure of the system, a very 
different temperature behavior is found. Specifically, in 
the new four-chain (case 3) regime we found a thermally 



induced structural phase transition, which shows up by 
the appearance of a plateau in the curves of the mean 
square displacement versus temperature. For the four- 
chain (case 2) regime an analysis of the pair-correlation 
function shows that the external chains, composed by 
particles with larger charge, exhibits a relatively ordered 
structure along the unconfined direction, even for a tem- 
perature larger than the melting temperature as deter- 
mined by the Lindemann-like criterion. The correlation 
function, indicates a qualitatively distinct melting in the 
X and y directions and for the distinct types of particles. 

To summarize, the two- to four-chain transition was 
found previously to be a first order transition (at T = 0)) 
and we found here that the density at which the transi- 
tion occurs is independent of temperature. The different 
melting temperatures exhibit small jumps at this phase 
transition. This is different for the different four-chain 
regimes, which was previously found to crystallize only 
in two configurations with a first order transition between 
them. Here we predict that a new four-chain phase exist 
intermediate between the previous two four-chain config- 
urations that we studied in Rcf. [10]. As a consequence, 
the transitions between the different four-chain phases 
are now continuous, i.e. second order, instead of first 
order transitions. In addition, the density at which the 
transition takes place depends on temperature which is 
a qualitative difference from the mono-disperse case 

For a « 1, the system can be found in a mixed config- 
uration which we studied in the regime of three chains. 
The configuration is characterized by a random distribu- 
tion of the two types of charges over the chains. In this 
regime, we found that an increase of temperature results 
in a more symmetric distribution of particles over the dif- 
ferent chains, and an increase of the melting temperature 
with density, except for high densities. 

The predicted melting phase diagram for the binary 
guasz-one-dimensional system can be investigated exper- 
imentally in systems as e.g. dusty plasmas and colloids. 
Previously, harmonic confinement potentials have been 
realized in dusty plasma p3 | and in colloidal [37| sys- 
tems. The phase diagram can be scanned experimentally 
by changing e.g. the density. Furthermore, in a dusty 
plasma it is also possible to vary the effective tempera- 
ture which is not possible for a colloidal system where the 
particles move in a liquid environment. But for the latter 
system one can also use paramagnetic particles which, in 
the presence of a perpendicular magnetic field, become 
magnetized resulting in a dipole inter-particle interaction 
potential. The strength of this potential can be tuned by 
the external magnetic field 25] . Because the obtained re- 
sults are generic, and qualitatively do not depend on the 
specific functional form of the repulsive interaction po- 
tential, similar melting behaviors are predicted for such 
a system of dipoles. Because in classical melting it is the 
thermal energy versus the average interaction energy that 
is the relevant parameter the temperature scale in Fig. 
3 can also be replaced by the inverse of the interaction 
energy. 
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Appendix A 



The zero-temperature phase diagram for a = 2 is 
shown in Fig. I12[ where we included the new four-chain 
(case 3) regime. A sketch of the four-chain (case 3) con- 
figuration is shown as inset in Fig. |6]and the parameters 
characterizing this configuration are plotted in Fig. [5] for 
K = 1 and a — 2. The phase transitions between the 
different four-chain configurations is of second order. 
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